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METHOD AND APPARATUS FOR GRADIENT 
ELECTROMAGNETIC INDUCTION WELL LOGGING 
BACKGROUND OF THE INVENTION 

The present invention relates in general to electromagnetic (EM) geo- 
physical method of detenriining the electrical conductivity (or resistivity) of 
earth formations penetrated by a well bore. The method can be applied 
for studying the underground geological structures in mineral, hydrocarbons 
and groundwater exploration, for example, for deterrnining the low conduc- 
tive zones associated with oil and gas reservoirs. 

Formation conductivity (or resistivity) determination from a wellbore is 
one of the oldest geophysical techniques to be applied in oil and gas explo- 
ration and production. The idea and principles of a geophysical method of 
EM induction well logging were introduced in the pioneering paper by H. 
G. Doll, Introduction to Induction Logging of Wells Drilled with Oil Based 
Mud, Journal of Petroleum Technology, vol. 1, p. 148, Society of Petroleum 
Engineers, Richardson lex. (1949). Over the decades, many modifications 
and improvements have been made in this area. These modifications can be. 
found in the numerous patents, for example, U.S. Pat. No. 5.452.762 issued 
to Beard et al., U.S. Pat. No. 5,781,436 issued to Forgang et al., U.S. Pat. 
No. 6,556,016 B2 issued to Gao et al. 

In conventional induction well logging, an instrument having transmitter 
coils and receiver coils substantially parallel to the borehole axis is lowered 
into the borehole and measures the magnetic field generated by the eddy cur- 
rents induced in the earth formations. Certain earth formations containing 
petroleum and permitting the petroleum to flow through the rock comprising 
the formation have certain physical characteristics well known in the art. For 
example, hydrocarbon-bearing formations are typically very resistive, while 
water-bearing formations are typically very conductive. A limitation to the 
EM induction well logging method known in the art is that the response 
of the typical EM induction logging instrument is largely dependent on the 
conductive layers, while the effect of the nonconductive hydrocarbon-bearing 
reservoirs can be masked by the influence of the conductive layers. 

At the same time, it is well known that the rock formations surrounding 
the borehole may be anisotropic with regard to the conduction of electrical 
currents. The phenomenon of electrical anisotropy is generally a consequence 
of either microscopic or macroscopic geometry, or a combination thereof. The 



reservoir anisotropy was identified as the major factor affecting the induction 
well logging data by Klein et al. (1997). Basically, the authors of the cited 
paper found that in a hydrocarbon-bearing reservoir at least two separate 
components of resistivity can influence the induction instrument: the resis- 
tivity measured with current flowing parallel to the bedding planes, which 
is called the transverse or horizontal resistivity pn, and whose inverse is the 
horizontal conductivity, <Th = l/p*; ^ e resistivity measured with a current 
flowing perpendicular to the bedding plane, which is called the longitudi- 
nal or vertical resistivity, p Vi a nd who se inverse is the vertical conductivity 
a v = l/p v . The ration, A = \Jcrh/vv> is called the anisotropy coefficient of 
the medium. Klein et al. (1997) found that the anisotropy coefficient for 
a hydrocarbon-bearing reservoir could rich the value up to 100. In this sit- 
uation the conventional induction well logging tool, oriented substantially 
perpendicular to the bedding planes, is sensitive to the horizontal resistivity 
of the formation only, while the instrument oriented at an angle with the 
bedding planes reads the apparent conductivity response, a a , which can be 
any value between y/a^cTv and ov 

One solution to this limitation known in the art has emerged in recent 
years. It is based on introducing transverse transmitter coils and receiver 
coils, with magnetic moments oriented perpendicular to the borehole axis. 
The idea of transverse induction coil measurements appeared first in the for- 
mer Soviet Union (Eidman, 1970; Kaufman and Kaganskii, 1972; Tabarovsky 
et al., 1976). It has received further development in many inventions. For- 
gang et al. (U.S. Pat. No. 5,781,436) introduced method and appara- 
tus for transverse electromagnetic induction well logging. Wu (U.S. Pat. 
No. 5,886,526) described a method of determining anisotropic properties 
of anisotropic earth formations using a multi-spacing induction tool. Hagi- 
wara (U.S. Pat. No. 5,966,013) disclosed a method of determining certain 
anisotropic properties of formation. Gao et al. (U.S. Pat. No. 6,393,364) 
introduced an iterative method for determining the horizontal and vertical 
resistivity. Gupta et al. (U.S. Pat. No. 5,999,883) considered a triad in- 
duction tool. Gao et al. (U.S. Pat. No. 6,556,016 B2) introduced an induc- 
tion method for determining the dip angle of an anisotropic earth formation 
surrounding a wellbore. In yet another development Zhdanov, Kennedy, 
and Peksen, in the article entitled "Foundations of the tensor induction well 
logging," developed the basic physical and mathematical principles of elec- 
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tromagnetic tensor induction well logging in anisotropic formation. 

Another limitation of the induction well logging data, known in the art, is 
the significant borehole effect that can distort the response of the formation 
surrounding the borehole. This effect has to be eliminated from the data, 
otherwise interpretation would be erroneous. The methods for the reduc- 
tion of the wellbore effect were developed by Tabarovsky and Epov (1972). 
They introduced a "frequency focusing" technique based on induction volt- 
age measurements at two frequencies and combining the observed signal in 
a manner so that the effects of eddy currents flowing within the wellbore 
can be substantially eliminated from the final result. This "dual frequency" 
signal is widely used in interpretation of conventional induction logging data 
and in transverse induction coil measurements as well (Forgang et al., U.S. 
Pat. No. 5,781,436). 

The above references are incorporated herein by reference. 

The foregoing attempts to determine vertical and horizontal resistivity 
around a deviated borehole and/or the dip angle of the formation met with 
varying degrees of success. However, there remains a need for improved 
methods and apparatus for high resolution hydrocarbon-bearing reservoir 
identification and characterization. A new technique is therefore needed. 

SUMMARY OF THE INVENTION 

The present invention provides a new method and apparatus for deter- 
mining the conductivity of anisotropic formations surrounding a borehole. 
The method comprises measuring the gradient of the magnetic field by an 
electromagnetic gradient induction logging tool in a borehole in the subsur- 
face formations. The method of the invention is capable of measuring the 
conductivity distribution of earth formations penetrated by a wellbore with 
higher resolution than conventional induction logging data, because the gra- 
dient of the magnetic field is significantly more sensitive to the conductivity 
variations than the magnetic field itself. 

An embodiment of the invention provides an apparatus which includes 
at least one transmitter coil and two closely positioned parallel receiver coils 
separated by a small vector Ar, with the magnetic moment direction of the 
transmitter parallel or different from the magnetic moment direction of the 
receivers. The transmitter generates a harmonic (frequency domain) primary 
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electromagnetic (EM) field which propagates through the anisotropic forma- 
tion surrounding the borehole. The EM field scattered by the formations 
is recorded by the receivers. The pair of closely located receivers with the 
parallel magnetic moments measures the gradient in the direction Ar of the 
component of the magnetic field, parallel to the receiver moment direction. 
A central processing unit collects the recorded amplitudes and phases of the 
EM field scattered by the formation and determines the EM parameters of 
the formation, such as horizontal and vertical conductivities and the dip 
angle of the bedding. 

In another embodiment of the invention, the frequency domain current in 
the transmitter is generated for at least several frequencies, and the receivers 
measure the signal at several frequencies. The multi-frequency gradient mea- 
surements are used for frequency gradient EM sounding of the medium at 
different distances from the wellbore. In addition, the multi-frequency mea- 
surements aTe combined in the gradient magnetic field signal with a reduced 
borehole effect. 

Another embodiment of the invention provides a new method and appa- 
ratus for measuring the conductivity of anisotropic formations surrounding a 
borehole. The apparatus comprises at least one receiver coil and two closely 
positioned parallel transmitter coils, separated by a small vector Ar, with the 
magnetic moment direction of the transmitters parallel or different from the 
magnetic moment direction of the receiver. The two transmitters generate 
successively a harmonic (frequency domain) primary EM field which prop- 
agates through the anisotropic formation, surrounding the borehole. The 
EM field scattered by the formations is recorded by the receiver. Due to 
the reciprocity principal, the difference of the induced voltages in the re- 
ceivers, measured successively for the signals from the two transmitters, is 
proportional to the gradient in the direction Ar of the component of the 
magnetic field, parallel to the receiver moment direction. A central process- 
ing unit collects the recorded amplitudes and phases of EM field scattered 
by the formation and determines the EM parameters of the formation, such 
as horizontal and vertical conductivities and the dip angle of the bedding. 

In another embodiment, the frequency domain current in the transmit- 
ters is generated at several frequencies. The multi-frequency gradient mea- 
surements are used for frequency gradient EM sounding of the medium at 
different distances from the wellbore. In addition, the multi-frequency mea- 
surements are combined in the gradient magnetic field signal with a reduced 
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borehole effect. 

Another embodiment of the invention provides a new method and appa- 
ratus for measuring the conductivity of an anisotropic formation surrounding 
a borehole. The apparatus comprises a tri-axial EM induction transmitter 
and a pair of closely located tri-axial EM induction receivers, separated by 
a small vector, Ar. One set of tri-axial receivers detects three components 
of the magnetic field due to each of three transmitters for a total of nine 
signals, which can be mathematically organized as a tensor array which we 
call the induction tensor. The entire instrument comprising a pair of tri-axial 
receivers detects a gradient in the direction Ar of each of three components 
of the magnetic field in the receivers due to each of the three transmitters 
for a total of nine gradient signals. A central processing unit collects the 
recorded amplitudes and phases of the EM field scattered by the formation 
and determines the EM parameters of the formation, such as horizontal and 
vertical conductivities and the dip angle of the bedding. 

Another embodiment of the invention provides a new method and appa- 
ratus for measuring the conductivity of an anisotropic formation surrounding 
a borehole. The apparatus comprises a pair of closely located tri-axial EM 
induction transmitters, separated by a small vector, Ar, and a tri-axial EM 
induction receiver. Due to the reciprocity principle, this instrument detects 
a gradient in the direction Ar of each of three components of the magnetic 
field in the receiver due to each of the three transmitters for a total of nine 
gradient signals. A central processing unit collects the recorded amplitudes 
and phases of the EM field scattered by the formation and determines the EM 
parameters of the formation, such as horizontal and vertical conductivities 
and the dip angle of the bedding. 

In another embodiment the frequency domain current in the transmitters 
is generated at several frequencies. The multi-frequency measurements in 
the tri-axial receivers are combined in the nine gradient magnetic field sig- 
nals which are used for the frequency gradient EM sounding of the medium at 
different distances from the weilbore. In addition, the multi-frequency mea- 
surements are combined in the gradient magnetic field signal with a reduced 
borehole effect. 

Another embodiment of the invention provides a new method and appa- 
ratus for measuring the conductivity of an anisotropic formation surrounding 
a borehole. The apparatus comprises a tri-axial EM induction transmitter 
and nine pairs of closely located parallel receiver coils. Three pairs of the 
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parallel receivers have the same orientation along the axis x, while they are 
separated by small distances in the directions x, y, and z, respectively, where 
the x, y, and z axes correspond to the coordinate system of the tool. An- 
other three pairs of the parallel receivers have the same orientation along the 
axis y, and they are separated by small distances in the directions x, y, and 
z respectively. The last three pairs of the parallel receivers have the same 
orientation along the axis z, and they are separated by small distances in 
the directions x, y, and z respectively. This instrument detects a gradient of 
the induction tensor, which is a tensor of rank 3 with a total number of 27 
gradient components. A central processing unit would collect the recorded 
amplitudes and phases of the EM field scattered by the formation and de- 
termine the EM parameters of the formation, such as horizontal and vertical 
conductivities and the dip angle of the bedding. 

In another embodiment the frequency domain current in the transmitters 
is generated at several frequencies. The multi-frequency measurements in 
the tri- axial receivers are combined in the 27 gradient magnetic field signals 
with a reduced borehole effect. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1(a) illustrates the coil configuration of a gradient induction tool 
consisting of one transmitter coil and two closely positioned parallel receiver 
coils separated by a small vector Ar. 

FIG. 1(b) shows the coil configuration of a reciprocal gradient induction 
tool consisting of one receiver coil and two closely positioned parallel trans- 
mitter coils separated by a small vector Ar. 

FIG. 2(a) presents the coil configuration of a tri-axial gradient induction 
tool consisting of a tri-axial EM induction transmitter and a pair of closely 
located tri-axial EM induction receivers, separated by a small vector Ar. 

FIG. 2(b) shows the coil configuration of a reciprocal tri-axial gradient 
induction tool consisting of a tri-axial EM induction receiver and a pair of 
closely located tri-axial EM induction transmitters, separated by a small 
vector Ar. 

FIG. 3 presents a schematic diagram of an instrument with an arrayed 
induction coil configuration consisting of a tri-axial transmitter array and 
three pairs of receiver coils; the pair of receiver coils E} z and R* is measuring 
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the dH%/dx (/? = x,y,z) components; the pair of receiver coils B^. and B% is 
measuring the dH^/dz, (p = z) components; the pair of receiver coils 

and E* is measuring thedHfj/dx, (J3 = x,y, z) components. 

FIG. 4 presents a schematic diagram of an instrument with an arrayed 
induction coil configuration consisting of three mutually orthogonal trans- 
mitter coils, displaced along the z-axis and oriented in the x- t y- } and jar- 
directions, and three pairs of receiver coils; the pair -of receiver coils R\ and 
R% is measuring the dH%/dx (/3 = z,y, z) components; the pair of receiver 
coils R]. and E* is measuring the dH£/dz> (0 = x, y, z) components; the pair 
of receiver coils -Rj and E% is measuring the dH^/dx, (P = x t y,z) compo- 
nents. 

FIG. 5 presents a schematic diagram of an instrument with an arrayed 
induction coil configuration consisting of three mutually orthogonal trans- 
mitter coils, displaced along the z-axis and oriented in the y-, and z- 
directions, and three pairs of receiver coils; the pair of receiver coils R\ and 
E% is measuring the dH%/dz (0 = z) components; the pair of receiver 
coils R* and R% is measuring the dHj*/dz, {fi — x, y, z) components; the pair 
of receiver coils R^ and i?£ is measuring the dH^/dz, (/? = x, y, z) compo- 
nents. 

FIG. 6 illustrates the tensor gradient induction instrument in an earth 
formation with several dipping beds with different vertical and horizontal 
resistivities, where the a/, ]/, and z* axes correspond to the coordinate system 
of the instrument, and the x t y i and z axes correspond to the coordinate 
system of the anisotropic layered earth formation. 

FIG. 7 shows the graphs of the theoretical gradient induction logging 
curves for a three-layer anisotropic model with the middle layer with the 
horizontal resistivity = 1 Ohm-m, vertical resistivity p v = 100 Ohm-m, 
and the thickness of 0.1 m. The horizontal and vertical resistivities of the 
upper and lower layers are equal to 1 Ohm-m, and 5 Ohm-m, respectively. 
Panel (a) presents horizontal resistivity variation versus depth and panel 
(b) shows vertical resistivity variation versus depth, wherein the bold solid 
lines display the true resistivity profile and the thin solid line shows the 
corresponding apparent resistivity profiles. The dual frequency H" d and 
data are presented in panels (c) and (d), respectively; panels (e) and (f) 
illustrate the dual frequency gradient components versus depth, dH*Jdz 
and dH% d /dz, respectively. The operating frequencies are 20 and 200 kHz. 

FIG. 8 shows the graphs of the theoretical gradient induction logging 
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curves for a three-layer anisotropic model with the middle layer having a 
horizontal resistivity of ph = 10 Ohm-m, a vertical resistivity of p v = 100 
Ohm-m, and a thickness of 0.5 m. The horizontal and vertical resistivities of 
the upper and lower layers are equal to 1 Ohm-m and 5 Ohm-m, respectively. 
Panel (a) presents horizontal resistivity variation versus depth and panel 
(b) shows vertical resistivity variation versus depth, wherein the bold solid 
lines display the true resistivity profile and the thin solid line shows the 
corresponding apparent resistivity profiles. The dual frequency H' d and 
data are presented in panels (c) and (d), respectively; panels (e) and (f) 
illustrates the dual frequency gradient components versus depth, BH^Jdz 
and dH^/dz, respectively. The operating frequencies are 20 and 200 kHz. 

FIG. 9 shows the graphs of the theoretical gradient induction logging 
curves for a three-layer anisotropic model with the middle layer having hori- 
zontal resistivity of ph = 1 Ohm-m, a vertical resistivity of p v = 100 Ohm-m, 
and a thickness of 2 m. The horizontal and vertical resistivities of the upper 
and lower layers are equal to 1 Ohm-m and 5 Ohm-m, respectively. Panel (a) 
presents horizontal resistivity variation versus depth and panel (b) shows ver- 
tical resistivity variation versus depth, wherein the bold solid lines display 
the true resistivity profile and the thin solid line shows the corresponding 
apparent resistivity profiles. The dual frequency H% d and data are pre- 
sented in panels (c) and (d), respectively; panels (e) and (f) illustrate the 
dual frequency gradient components versus depth, dH' d /dz and dH^/dz^ 
respectively. The operating frequencies are 20 and 200 kHz. The observed 
and predicted data are shown by dotted and solid lines in panels (c) and (d), 
(e) and (f ) , respectively. 

FIG. 10 shows the graphs of the theoretical gradient induction logging 
curves for a three-layer anisotropic model with the middle layer having hori- 
zontal resistivity of ph = 10 Ohm-m, a vertical resistivity of p„ = 100 Ohm-m, 
and a thickness of 2 m. The horizontal and vertical resistivities of the upper 
and lower layers are equal to 1 Ohm-m and 5 Ohm-m, respectively. Panel (a) 
presents horizontal resistivity variation versus depth and panel (b) shows ver- 
tical resistivity variation versus depth, wherein the bold solid lines display 
the true resistivity profile and the thin solid line shows the corresponding 
apparent resistivity profiles. The dual frequency and data are pre- 
sented in panels (c) and (d), respectively; panels (e) and (f) illustrates the 
dual frequency gradient components versus depth, 3H*Jdz and dH£ d /dz, 
respectively. Operating frequencies are 20 and 200 kHz. The observed and 
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predicted data are shown by dotted and solid lines in panels (c) and (d), (e) 
and (f), respectively. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

The conceptual sketches for coil arrangements for downhole gradient in- 
duction instruments are shown in FIGS. 1, 2, 3, 4 & 5. 

FIG. 1(a) illustrates a preferred embodiment in which a gradient induc- 
tion instrument consists of one transmitter coil with the magnetic moment 71, 
and two closely positioned parallel receiver coils with the magnetic moments 
JR^ and B^ n) separated by a small vector, Ar, with the magnetic moment 
direction of the transmitter, T t , parallel or different from the magnetic mo- 
ment direction of the receivers. In the preferred embodiment, it is assumed 
that the distance between two parallel receivers is much smaller than the dis- 
tance L between the transmitter and a center of the vector Ar, connecting 
two receivers: |Ar| << L. This requirement ensures that the receiver pair 
measures the gradient of the magnetic field component H m at the direction 
of the vector Ar : 

dn^" |Ar| ' W 

where H m and are the components of the magnetic field, parallel to the 
receiver moment direction, which are measured by receivers and H^, 
respectively, and nAr is a unit vector in the direction of vector Ar. 

FIG. 1(b) illustrates a preferred embodiment in which a reciprocal gra- 
dient induction tool consists of one receiver coil with the magnetic moment 
Ri t and two closely positioned parallel transmitter coils with the magnetic 
moments T m and separated by a small vector Ar, with the magnetic 
moment direction of the receiver, Ri parallel or different from the magnetic 
moment direction of the transmitters. In the preferred embodiment, it is 
assumed that the distance between two parallel transmitters is much smaller 
than the distance L between the receiver and a center of the vector Ar, con- 
necting two transmitters: |Ar| « L. The transmitters T m and operate 
successively, so that the receiver Ri measures the magnetic fields, generated 
by each of these transmitters separately. Due to the reciprocity principle, the 
difference between the two measurements of the magnetic field in the receiver 
due to two different transmitters is equal to the gradient of the magnetic field 
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component H m at the direction of the vector Ar, measured by receivers 
and of the coil configuration shown in FIG. 1(a). 

FIG. 2(a) illustrates the preferred embodiment in which the instrument 
consists of a tri-axial EM induction transmitter and a pair of closely located 
tri-axial EM induction receivers, separated by a small vector Ar. In the 
transmitter array, there is a triad of transmitter coils, T Xi T y and T z , each 
oriented along a respective axis. In each receiver array a similarly oriented 
triad of receiver coils, R^ t #J, and R\ and R% y i^, and Rl, is provided. It 
is assumed that the three coils in each triad represent actual coils oriented 
in mutually perpendicular directions, with the z-axis corresponding to the 
axis of the borehole. However, this coil arrangement is not a requirement 
of the invention, as any coil orientation can be synthesized by performing a 
suitable rotation on coils with different triad orientations. Such rotations are 
described in depth in Zhdanov et al. (2001), which is incorporated herein by 
reference. 

In a preferred embodiment, each of the coils in the transmitter triads !T X , 
T yt and T g is parallel to the corresponding coils in the receiver triad R*, 
and R\, and/or i2£, and Rl and each is displaced from the center of the 
vector Ar, connecting two receiver triads by a distance L. 

One set of tri-axial receivers detects three components of the magnetic 
field due to each of three transmitters for a total of nine signals, which can 
be mathematically organized as a tensor array which we call the induction 
tensor, H. This tensor can be written in matrix notation as 



HI H% HI 

h; hi H$ 



where superscripts refer to the transmitter components and subscripts refer 
to the receiver components. In the preferred embodiment, it is assumed 
that the distance between two triads of receiver coils is much smaller than 
the distance L between the transmitter triad and a center of the vector Ar, 
connecting two receiver triads: |Ar| « L. This requirement ensures that 
the receiver triads measure the gradient of the magnetic field components at 
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the direction of the vector Ar, which form a gradient tensor 



r am am am n 



g = ^S = 28* j%£ 



(2) 



dn^ %%r f$r 




where 




(3) 



and and iff 2 (a,/? = x } y, z) are the components of the magnetic field, 
which are measured by receivers and R\ and i? 2 , i? 2 , and B%, respec- 

tively, and n^r is a unit vector in the direction of vector Ar. 

An alternative embodiment is shown in FIG. 2(b) such that the tool 
consists of a tri-axial EM induction receiver and a pair of closely located 
tri-axial EM induction transmitters, separated by a small vector Ar. In the 
receiver array, there is a triad of receiver coils R^, Ry, and i^, each oriented 
along a respective axis. In each transmitter array a similarly oriented triad 
of transmitter coils Ti, I*, and T% and 7%, 2* and 2* is provided. It 
is assumed that the three coils in each triad represent actual coils oriented 
in mutually perpendicular directions, with the z-axis corresponding to the 
axis of the borehole. However, this coil arrangement is not a requirement 
of the invention, as any coil orientation can be synthesized by performing 
a suitable rotation on coils with different triad orientations. Due to the 
reciprocity principle, the difference between the two measurements of the 
magnetic field in the receivers due to two different transmitters triads is 
equal to the gradient of the magnetic field components at the direction of 
the vector Ar, measured by receiver triads i£, i£ and R l zi and R 2 x1 R* and 
Rl of the coil configuration shown in FIG. 2(a). , 

While particular embodiments of the present invention have been dis- 
closed, it is to be understood that various different modifications are pos- 
sible and are contemplated within the scope of the specification, drawings, 
abstract and appended claims. FIGS 3, 4, and 5 provide conceptual sketches 
for several different coil arrangements for the gradient tensor induction well 
logging (GTIWL) downhole instruments, which are included in the preferred 
embodiments of this invention. 

FIG 3 illustrates a preferred embodiment in which the tool consists of 
a tri-axial transmitter array and three pairs of receiver coils. The pair of 
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receiver coils R\ and R* is measuring the dH%/dx « (iff 2 - iff 1 ) /Ax, 0 = 
a;, y, components. The pair of receiver coils i£ and R* is measuring the 
dH%/dz » (iff 2 —Hj*ty jAz % # = components; the pair of receiver 

coils and is measuring the dH$/dx « (iff - H&) /Arc, /? = z, 
components. 

FIG 4 illustrates an alternative embodiment in which the tool consists 
of three mutually orthogonal transmitter coils, displaced along the z-axis 
and oriented in the x-, y-, and z- directions, and three pairs of receiver 
coils. The pair of receiver coils R\ and R% is measuring the DH^/dx « 
(iff 2 — iff 1 ) /Ax, p — x,2/,z, components. The pair of receiver coils R* 

and it^ is measuring the dH%/dz « (iff 2 - Hg) /Az, 0 = x y y,z, com- 
ponents; the pair of receiver coils f?J and R* is measuring the dH^jdx « 

(Hfj 2 - Hg 1 ) /Ax, /? = x, y, 2, components. 

FIG 5 illustrates an alternative embodiment in which the tool consists 
of three mutually orthogonal transmitter coils, displaced along the 2:- axis 
and oriented in the x-, y-, and z- directions, and three pairs of receiver 
coils. The pair of receiver coils R\ and R* is measuring the dH%/dz « 
(iff 2 — iff 1 ) / Az 7 = x,y t 2, components. The pair of receiver coils R%. 

and R* is measuring the dHg/dz « (iff 2 - iff*) /Az y - x, y, z, com- 
ponents; the pair of receiver coils R^ and R* is measuring the dHfj/dz « 

(H$ 2 - ffjf 1 ) /Az, /3 = x,y t z, components. 

In another preferred embodiment the tool consists of a tri-axial transmit- 
ter array and up to 27 pairs of receiver coils, measuring some or all magnetic 
gradient components: dH%/&y (a,fi J j = x i y 1 z). In yet another preferred 
embodiment the tool consists of three mutually orthogonal transmitter coils, 
displaced along the z-axis and oriented in the x-, y-, and z- directions, and 
up to 27 pairs of receiver coils, measuring some or all magnetic gradient 
components: dH%/dy (or, /?, 7 = x, y ) z) . 

We should notice that it is weD known in the art of exploration geophysics 
that the gradient of the magnetic field is significantly more sensitive to the 
local variations of the conductivity than the magnetic field itself. Thus, the 
gradient tensor induction well logging (GTIWL) measurements are capable 
of measuring the conductivity distribution of earth formations penetrated by 
a wellbore with a higher resolution than conventional induction logging data. 
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For interpretation of the instrument measurements, we, generally, con- 
sider an appropriate formation model. The interpretation problem is formu- 
lated for the GTIWL data collected in a deviated borehole within a planar 
stratified medium with electrical anisotropy in each layer. We assume that a 
model of the horizontally layered medium is characterized in the coordinate 
system {rr, y, z} by 1-D piecewise-constant distribution of the conductivity 
along the z direction. Each layer has transverse isotropic (TT) conductivity 
defined by the conductivity tensor ai(z) 




0 0 

< o 

0 K 

and the corresponding piecewise-constant anisotropy coefficient 

* = i=l,.,N. 
The positions of the layer's boundaries are given by equations 



(4) 



(5) 



z = Zi, t=l,..,JV*-l. ! (6) 

We consider first the response of a tri-axial gradient induction instrument 
shown in FIG 2(a) in an unbounded, homogeneous, transversely isotropic 
(TI) conductive medium. Every receiver triad of this instrument detects 
three components of the magnetic field due to each of three transmitters for 
a total of nine signals: 



•#x(i) 

■^vft) ^tt) H v(i) 

H z(l) #J(l) J 



and 



H< 2 > = 



TJX 
TJX 

#z(2) 



^(2) 
ti z(2) 



^x(2) 

H y{2) 
H%(2) J 



where superscripts indicate the transmitter components and subscripts rep- 
resent the receiver components, and index (1) or (2) corresponds to the first 
or the second receiver triad. Using formula (3), these measurements can be 
transformed into the gradient tensor G described by formula (2). 
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In the Cartesian system of coordinates (x 1 y i z) i with the axis z directed 
along the axis of symmetry of the TI conductive medium, the expressions for 
the induction tensor components are written as (Zhdanov et al., 2001) 

e^* \iuiJLQ<Jh ikh3 — k h kvX 2 21^7? 



TTX 



As 



* J 



47T 



ik h T - k 2 x 2 2ik h x 2 ik h 



--f+ 



• (A&g 2 -f 1) Zik h x 2 Zx 2 ] 
r* r 4 r^J' 



Ht=m = xy 
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,ik h r 
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(7) 



(8) 



(9) 



(10) 



(11) 



H ' ■ p** + 7" — ? ^- + h-J' < 12 ) 

where we use the following notations: p = + 2/ 2 > 5 = vV 2 + A'^ 2 , 
A 2 = <Jh/cr V) r = vV 2 + — ivfXTh) K = iujfia v , u) is the circular 

frequency of the current in the transmitter, cr* is the horizontal component 
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of the conductivity and a v is the vertical component of the conductivity. 
By differentiating the expressions (7)-(12) in the direction n^r* we find the 
components of the gradient tensor G. 

The magnetic field components are given in formulae (7)- (12) and, cor- 
respondingly, in expression (2) for the gradient tensor in a coordinate sys- 
tem defined by the horizontal and vertical principal axes of the transverse 
isotropic media. The oil and gas industry widely uses deviated wells and di- 
rectional drilling. That is why it is important to study the GTIWL response 
in a deviated well, where the orientation of the transmitter and receiver coils 
will be arbitrary with respect to this coordinate system. In order to use the 
representation of the gradient tensor G for an instrument located in an ar- 
bitrary orientation with respect to the tensor principal axes, it is necessary 
to transform the transmitter moment in the instrument frame (denoted by 
(x^j/,^)) into the medium coordinates (denoted (x,y, z)), as it is shown in 
FIG 6. This transformation can be made by application of the rotational 
matrix R, given by 

cos a cos P cos a sin ft —sin a 

— sin/? cos/3 0 
sin a cos /? sinasin/? cos or 

With the use of the rotational matrix, the relations between the matrix 
G' of the gradient tensor in the instrument coordinates and its matrix G 
in the coordinates defined by the conductivity tensor principal axes can be 
described by the following formula: 



R = 



(13) 



G' = RGR- 1 = BjSSF. (14) 

We should note, however, that both the conventional induction logging 
data and GTIWL data are strongly affected by the eddy currents flowing 
in the wellbore. This effect must be eliminated from the data, otherwise 
interpretation would be erroneous. The methods for the reduction of the 
wellbore effect were developed in the pioneering paper by labarovsky and 
Bpov (1972). They introduced a "frequency focusing" technique based on 
induction voltage measurements at more than one frequency and combining 
the observed signal in a manner so that the effects of eddy currents flowing 
within the wellbore can be substantially eliminated from the final result. 
The further modifications of this technique can be found in U.S. Pat. No. 
5,781,436 issued to Pbrgang et al. 
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In order to remove the borehole and invasion zone effects from the gradient 
induction logging data, one can calculate the dual frequency signal for every 
component of the gradient tensor: 

G i {f) = G h (h)-(f£\G h {h), (15) 

where Gd (/) is a dual component (A/m), and fi and fa are operating fre- 
quencies (Hz). and G/ 2 (/2) are the corresponding components of 
the gradient tensor G given by equation (2), observed or calculated at these 
frequencies. As we can see a dual component can be calculated by measur- 
ing gradient of the magnetic field components at two various frequencies. In 
practice, the quadrature (imaginary) components are usually used for con- 
ductivity calculation. 

We will present below, as an example, the dual frequency signals for the 
different components of the gradient tensor (2) in a simple case of a vertical 
borehole within a horizontally layered anisotropic formation. In this case, 
the vertical magnetic field component (coaxial) can be given by the formula 
(Zhdanov et al., 2001): 

Hl = ^(l-ik h L), . (16) 

where kh = V^W^i k a horizontal wave number, (= 1/ph) (S/m) is a 
horizontal conductivity (or resistivity), \i = /io = 47rl0~ 7 (H/m) is the free- 
space magnetic permeability, a; is an angular frequency with u = 2wf , and 
L is the distance between transmitter and receiver. The coaxial magnetic 
dual component is written as: 

HU (fi, h) = (1 - " (gj ^£3 (1 - ■ (17) 

where = V^iA*o^ j kh 2 ~ v^2AW^> w i — 27r/i, and 0J2 — 2?r/2 . 

In a similar way we can write a horizontal (coplanar) component of mag- 
netic field in an anisotropic medium: 

e^ L I 14- A 2 o «\ 

^--M 1- ^"^^)' (18) 
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where A = \f^h/^v is an anisotropy coefficient, and a v is the vertical con- 
ductivity of the medium. The coplanar magnetic dual component is written 



as: 



e ikh i L I 1 



±2l 

2A 2 



) 



(19) 



1>T kh * L 



Using the same approach, one can write the corresponding dual signals 
for every component of the induction tensor. In the case of the deviated bore- 
hole and arbitrary orientation of the instrument with respect to the medium 
coordinate frame (with respect to the layered formations), the expressions for 
the magnetic field components can be obtained from corresponding formulas 
developed in the medium coordinate system, by application of the rotational 
transformation (14) (Zhdanov et aL, 2001). 

In the case of the gradient instrument, we measure the spatial derivatives 
of the corresponding magnetic field components. For example, the vertical 
derivative of the vertical coaxial magnetic dual component is written as: 



Gld (/i./a) = ^ [-3 + 3ifc hl L + fc 2 ,!, 2 ] 



(20) 



The vertical derivative of the coplanar magnetic dual component is given by 
the expression: 



(fu h) = [-3 + 3iA Al £ + kil? - 



1 + ^ t 7 r 3 T 3 



-( 



•3 + Zik h2 L + 



1 + A 2 
2A 2 



■ (21) 



It is important to understand the basic difference between magnetic in- 
duction well logging data measured by conventional instruments and gradient 
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magnetic data measured by the new gradient induction instruments of this 
invention. FIG. 7 shows the graphs of the theoretical induction logs and 
gradient induction logs for a three-layer anisotropic model with the middle 
layer having a horizontal resistivity of ph = 1 Ohm-m, a vertical resistivity 
of fa — 100 Ohm-m, and a thickness of 0.1 m. The horizontal and verti- 
cal resistivities of the upper and lower layers are equal to 1 Ohm-m and 5 
Ohm-m, respectively. Therefore, in this model there is no variation in the 
horizontal resistivity (see panel a), while there is a thin intermediate layer 
with the anomalous vertical resistivity (see panel b). The dual frequency 
and H* d data are presented in panels (c) and (e), respectively; panels (d) and 
(f) illustrate the dual frequency gradient components versus depth, dH^Jdz 
and dHZj/dz, respectively. The operating frequencies are 20 and 200 kHz. 
One can see that the standard induction log (panel c) and the vertical gradi- 
ent of this log (panel d) cannot determine the presence of the middle layer, 
because they both are sensitive to the horizontal resistivity only. The 
component, and specially its vertical gradient dH^/dz, clearly indicate the 
presence of an intermediate layer with high vertical resistivity (panels e and 



We can observe a similar effect in a three- layer anisotropic model with 
the middle layer having a horizontal resistivity of ph = 10 Ohm-m, a vertical 
resistivity of p v = 100 Ohm-m, and a thickness of 0.5 m, presented in FIG 8. 
In this model the intermediate layer is different from the background section 
in both the horizontal and vertical resistivities. Therefore, we observe the 
anomalous behavior in all four induction logs, which are shown in panels (c) 
- (f) of FIG 8. However, the most visible anomaly is observed in the log of 
the vertical gradient dH^/dz (panel f). 

One can see from expressions (20) and (21) that the dual frequency signals 
are the functions of the horizontal, a hi and vertical, a v , conductivities. Our 
goal is to find these conductivities from the observed data. This is obviously 
a nonlinear inverse problem. This problem can be solved by applying the 
regularized Newton method as described by Zhdanov (2002): 



f). 



r n = B(<r„) -g, 

H„ = F£F n + v n l, 
<rn+i =<r„-H^ 1 l^, 



(22) 
(23) 

(24) 
(25) 
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where B is a nonlinear operator of the forward modeling described by the 
corresponding formulas for dual frequency signals (see, for example, formulae 
(20) and (21) for a homogeneous anisotropic medium). The suitable forward 
modeling method for the purpose of calculations of the predicted GTIWL 
data B(<r„) is disclosed in Zhdanov et al.„ 2001. Vector g is formed by 
the observed data (dual frequency gradient magnetic field signal), <r n is a 
vector of model parameters (horizontal and vertical conductivities, ah and 
fj v , and the relative deviation angle a of the borehole) on the nth iteration, 
<r apr is the a priori model, r is a residual vector of the difference between the 
predicted, "B(cr n ) } and observed data; F w is a Frechet derivative matrix; 
is a quasi Hessian matrix; and l£ n is the regularized direction of the Newton 
method on the nth iteration. 

The regularization parameter v is updated on each iteration according to 
a progression of numbers: 

Vk = ^og*; k = 0, 1, 2, n; q > 0. (26) 

The first iteration of the Newton method is run with v = 0. The initial value 
of the regularization parameter (vq) is determined after the first iteration 
(<r i ) as the ratio: 

^ JlB^-gll 2 

For any number (v h ) we calculate the misfit || B (<r„J — g \\*. The optimal 
value of the parameter v is the number i^o, for which we have satisfied the 
misfit condition: 

l|B(<0-g||*=S, (27) 

where 8 is the level of noise in the observed data. 

There are three unknown model parameters: the relative deviation angle 
a and the horizontal and vertical conductivities, an and <x v . Therefore, the 
Prechet derivative matrix can be calculated directly by taking derivatives of 
the corresponding expressions for the dual frequency signals with respect to 
unknown parameters. The matrix is: 



F = 



SSL 



gey 
SQL 



da 
da 

da 
dG** 

da J 



(28) 
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In the Newton algorithm, the calculations start with an initial guess, 
which usually is set equal to the a priori model, <ro = Vapr, and then updated 
on each iteration according to (25). The method is fast and converges after 
four-six iterations. 

In a general case, the Newton routine solves for three parameters at each 
GTTWL observation point. In a preferred embodiment of the invention, the 
Newton method described above is used to determine the apparent horizontal 
and vertical conductivities of the medium, and the relative deviation angle 
from the GTTWL data. 

FIG 9 shows an example of the apparent resistivity determination from 
the GITWL data using the Newton method for a three-layer anisotropic 
model with the middle layer having a horizontal resistivity of ph = 1 Ohm- 
m, a vertical resistivity of p v = 100 Ohm-m, and a thickness of 2 m. We can 
see that the apparent vertical resistivity (thin solid line in panel b) describes 
well the true resistivity of the model. The predicted and observed data are 
shown by dotted and solid lines in panels (c) and (d), (e) and (f) , respectively. 
One can see that the predicted data fit the observed data very well. 

FIG 10 shows a similar example for a three-layer anisotropic model with 
the middle layer having a horizontal resistivity of pn ~ 10 Ohm-m, a vertical 
resistivity of p v = 100 Ohm-m, and a thickness of 2 m. In this case, the 
apparent resistivities provide some indication about the location and the 
electrical resistivity of the intermediate layer, but do not describe the true 
resistivity profiles well enough. We observe the typical polarization horns 
in the apparent vertical resistivity curve (panel b). In this situation the 
apparent resistivities can be used as the starting models for the iterative 
nonlinear inversion for the layered model parameters. 

Another embodiment of the claimed subject matter exists when the gra- 
dient magnetic field measurements are used for determining the electrical 
conductivity of the multilayered anisotropic formations with the borehole 
arbitrarily oriented with respect to the bedding plane. In this case vector 
cr of the model parameter distributions consists of logarithms of horizontal 
and vertical conductivities (lno^ and In cr^, i=l, .., N), where sample layers 
forming a 1-D gepelectrical inverse model: 

cr = [m(7M,ma t ,i,m^,mcr V 2,m<7/ l 3,lna' 1 ,3....ma^,m , 



and iV is a number of sample layers. The induction logging inverse problem 
consists of finding a 1-D distribution of model parameters (a) that corre- 
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spond to an observed discrete set of gradient tensor induction data (d). The 
solution of this problem is nonunique and ill-posed. In the preferred embodi- 
ment of the invention we apply the regularization method, given by Tikhonov 
and Arsenin (1977) and Zhdanov (2002), which is reduced to solving the fol- 
lowing minimization problem: 

P°{<t) =<f>{<r)+vS{<r) = min, (29) 

where /^(cr) is the parametric functional, <t>(a) is a misfit functional, v is 
a regularization parameter, and S{cr) is a stabilizing functional. The misfit 
functional is chosen as the following form, 



</>{<r) = ||w d (Bo- - g 



(30) 



where Wj is a data weighting matrix. We follow the ideas of the general 
method of sharp boundary inversion ( Zhdanov, 2002), and use the minimum 
gradient support stabilizing functional: 



S(cr) = S M gs{ct) = 



|grada| 



(|gradof + e 2 ) 



1/2 



(31) 



The solution of the rninimization problem (29) can be obtained by the same 
Newton method (22) - (25) outlined above. 

Another embodiment comprises the multifrequency measurements of the 
GTIWL data for conducting a frequency gradient EM sounding of the medium 
at different distances from the wellbore. In this case a selected volume of the 
medium surrounding the borehole is descretized in the large number of small 
cells having constant horizontal and vertical conductivities.. Vector or of 
the model parameter distributions consists of logarithms of horizontal and 
vertical conductivities (ln<7w and In a^, i=l,.., JV) of the different cells: 

where TV is a number of cells. The induction logging inverse problem consists 
of finding a 3-D distribution of model parameters (cr) that corresponds to an 
observed discrete set of the multi-frequency gradient tensor induction data 
(d). The solution of this problem is based on the same principles of the 

regularization theory outlined in the previous embodiment of this invention. 
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What is claimed is: 

1. An apparatus for measuring the conductivity of the anisotropic earth 
formations penetrated by a wellbore, comprising 

a) at least one transmitter coil and two closely positioned parallel receiver 
coils separated by a small vector, Ar, with the magnetic moment direction 
of the transmitter parallel or different from the magnetic moment direction 
of the receivers, where the distance between two parallel receivers is much 
smaller than the distance L between the transmitter and the center of the 
vector Ar, connecting two receivers: |Ar| « L; 

b) at least one receiver coil and two closely positioned parallel transmitter 
coils, separated by a small vector, Ar, with the magnetic moment direction 
of the transmitters parallel or different from the magnetic moment direction 
of the receiver, where the distance between two parallel transmitters is much 
smaller than the distance L between the receiver and the center of the vector 
Ar, connecting two transmitters : |Ar| << L\ 

c) a tri-axial EM induction transmitter formed by three mutually orthog- 
onal transmitter coils, and a pair of closely located tri-axial EM induction 
receivers, each formed by three mutually orthogonal receiver coils, separated 
by a small vector, Ar, wherein the distance between two tri-axial EM in- 
duction receivers is much smaller than the distance L between the tri-axial 
transmitter and the center of the vector Ar: | Ar| « L; 

d) a pair of closely located tri-axial EM induction transmitters, separated 
by a small vector, Ar, and a tri-axial EM induction receiver, wherein the 
distance between two tri-axial EM induction transmitters is much smaller 
than the distance L between the tri-axial receiver and the center of the vector 
Ar: |Ar| « L\ 

e) means for receiving voltages induced in said receiver coils; 

f ) means for measuring a difference between the said voltages in the dif- 
ferent pairs of the receivers, or due to the different pairs of the transmitters. 

2. The apparatus as defined in claim 1 further comprising: 

a) a tri-axial transmitter array and three pairs of receiver coils, measuring 
the dH^/dx, DH^/dz, and dH$/dx (J3 = x,y,z) components; 

b) a tri-axial transmitter array and three pairs of receiver coils, measuring 
the H%/dz t dHg/dz t and dH?/dz {/3 = x,y, z) components; 
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c) a tri-axial trans mitt ex array and up to 27 pairs of receiver coils, mea- 
suring some or all components: dH^/fry (a,/?, 7 = x t y t z) ; 

d) three mutually orthogonal transmitter coils, displaced along the z-aoris 
and oriented in the a?-, y-, and z- directions, and three pairs of receiver coils, 
measuring the dH%fdx, dH£/dz, and dH$/dx (P — x>y 7 z) components; 

e) three mutually orthogonal transmitter coils, displaced along the z-axis 
and oriented in the a;-, y-, and z- directions, and three pairs of receiver coils, 
measuring the dH£/dz, dH^jdz^ and dH^/dz {0 = x t y t z) components; 

f) three mutually orthogonal transmitter coila, displaced along the 2-axis 
and oriented in the x-,. y-, and z- directions, and up to 27 pairs of receiver 
coils, measuring some or all components: dH%/&y (a, /?, 7 = x> y, z) . 

3. A method for measuring the conductivity of the anisotropic earth 
formations penetrated by a wellbore, comprising: 

a) measuring a gradient of magnetic field between two closely positioned 
parallel receiver coils, wherein the harmonic (frequency domain) electromag- 
netic field is generated by at least one transmitter coil, and the receiver coils 
are separated by a small vector, Ar, with the magnetic moment direction of 
the transmitter coil parallel or different from the magnetic moment direction 
of the receiver coils, wherein the distance between two parallel receiver coils 
being much smaller than the distance L between the transmitter and the 
center of the vector Ar, connecting two receivers, |Ar| « L f to ensure that 
the receiver pair measures the gradient of the induction magnetic field; 

b) using said gradient electromagnetic induction logging instrument for 
obtaining a plurality of measurements for different positions of the instrument 
along the borehole; 

c) obtaining from the said plurality of measurements the horizontal and 
vertical conductivities of the medium, and the relative dip angle of the forma- 
tion by inverting the gradient electromagnetic induction data using a model 
of a layered anisotropic formation and regularization methods of inverse prob- 
lem solution. 

4. The method of claim 3, wherein the frequency domain current in the 
transmitter is generated for at least several frequencies, and the receivers 
measure the signal at several frequencies, wherein the multi-frequency gra- 
dient measurements are used for frequency gradient EM sounding of the 
medium at different distances from the wellbore to produce a volume image 
of the anisotropic conductivity distribution around the borehole. 

5. A method for measuring the conductivity of the anisotropic earth 
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formations penetrated by a wellbore, wherein the method comprises: 

a) measuring a gradient of magnetic field at the receiver coil position, 
wherein the harmonic (frequency domain) or pulse (time domain) electro- 
magnetic field is generated by two closely positioned parallel transmitter 
coils, separated by a small vector, Ar, with the magnetic moment direction of 
the transmitters parallel or different from the magnetic moment direction of 
the receiver, wherein two transmitters generate successively a harmonic (fre- 
quency domain) primary EM field which propagates through the anisotropic 
formation, surrounding the borehole. 

b) using said gradient electromagnetic induction logging instrument for 
obtaining a plurality of measurements for different positions of the instrument 
along the borehole; 

c) obtaining from the said plurality of measurements the horizontal and 
vertical conductivities of the medium and the relative dip angle of the forma- 
tion by inverting the gradient electromagnetic induction data, using a model 
of layered anisotropic formation and regularization methods of the inverse 
problem solution. 

6. The method of claim 5, wherein the frequency domain current in the 
transmitter is generated for at least several frequencies, and the receivers 
measure the signal at several frequencies, wherein the multi-frequency gra- 
dient measurements are used for frequency gradient EM sounding of the 
medium at different distances from the wellbore to produce a volume image 
of the anisotropic conductivity distribution around the* borehole. 

7. A method for measuring the conductivity of the anisotropic earth 
formations penetrated by a wellbore, wherein the method comprises: 

a) conveying an electromagnetic logging instrument into a borehole, said 
logging instrument including a tri-axial EM induction transmitter and a pair 
of closely located tri-axial EM induction receivers, separated by a small vec- 
tor, Ar, wherein one set of tri-axial receivers detects three components of 
the magnetic field due to each of three transmitters for a total of nine signals 
forming an induction tensor; 

b) said logging instrument measuring a gradient in the direction Ar of 
each of three components of the magnetic field in the receivers due to each 
of three transmitters for a total of nine gradient signals forming a gradient 
tensor. 

c) using said gradient tensor induction well logging (GTTWL) instrument 
for obtaining a plurality of measurements for the different positions of the 
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instrument along the borehole; 

d) obtaining from the said plurality of measurements the horizontal and 
vertical conductivities of the medium and the relative dip angle of the forma- 
tion by inverting the gradient tensor induction well logging (GTIWL) data 
using a model of layered anisotropic formation and regularization methods 
for the inverse problem solution. 

8. The method of claim 7, wherein the frequency domain current in the 
tri-axial transmitter is generated for at least several frequencies, and the tri- 
atrial receivers measure the signal at several frequencies, wherein the multi- 
frequency measurements are combined in the gradient tensor induction signal 
to be used for frequency gradient EM sounding of the medium at different 
distances from the wellbore to produce a volume image of the anisotropic 
conductivity distribution around the borehole. 

9. A method for measuring the conductivity of the anisotropic earth 
formations penetrated by a wellbore, wherein the method comprises: 

a) conveying an electromagnetic logging instrument into a borehole, said 
logging instrument including a pair of closely located tri-axial EM induction 
transmitters, separated by a small vector, Ar, and a tri-axial EM induction 
receiver, wherein a set of tri-axial receivers detects three components of the 
magnetic field due to each of three transmitters for a total of nine signals 
forming an induction tensor; 

b) measuring by the said logging instrument a gradient in the direction 
Ar of each of three components of the magnetic field in the tri-axial receiver 
for a total of nine gradient signals forming a gradient tensor. 

c) using said gradient tensor induction well logging (GTIWL) instrument 
for obtaining a plurality of measurements for the different positions of the 
instrument along the borehole; 

d) obtaining from the said plurality of measurements the horizontal and 
vertical conductivities of the medium and the relative dip angle of the forma- 
tion by inverting the gradient tensor induction well logging (GTIWL) data 
using a model of layered anisotropic formation and regularization methods 
for the inverse problem solution. 

10. The method of claim 9, wherein the frequency domain current in the 
tri-axial transmitters is generated for at least several frequencies, and the tri- 
axial receivers measure the signal at several frequencies, wherein the multi- 
frequency measurements are combined in the gradient tensor induction signal 
to be used for frequency gradient EM sounding of the medium at different 
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distances from the wellbore to produce a volume image of the anisotropic 
conductivity distribution around the borehole. 



